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Abstract 

In this talk I give an elementary introduction to the key algorithm used in recent ap- 
plications of computational algebraic geometry to the subject of string phenomenology. I 
begin with a simple description of the algorithm itself and then give 3 examples of its use in 
physics. I describe how it can be used to obtain constraints on flux parameters, how it can 
simplify the equations describing vacua in 4d string models and lastly how it can be used to 
compute the vacuum space of the electroweak sector of the MSSM. 
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1 Introduction 



There is currently a great deal of interest in applying the methods of computational algebraic geometry 
to string phenomenology and closely related sub- fields of theoretical physics. For some examples of 
recent work see 121 131 IH IS1 IB1 13 IB1 ISl UDl UH U21 H31 and references therein. These papers utilise 
advances in algorithmic techniques in commutative algebra to study a wide range of subjects including 
various aspects of globally supersymmetric gauge theory [B El EJ [H O [B] , finding flux vacua in string 
phenomenology (7J [HJ [9], studying heterotic model building on smooth Calabi-Yau in non-standard 
embeddings [TO] and more besides [TT | \12 \ \13 \ [13]. 

Despite the wide range of physical problems which have been addressed within this context, the 
computational tools which are being used are all based, finally, on the same algorithm. The Buchberger 
algorithm [15] is at once what lends these methods their power and also the rate limiting step - placing 
bounds on the size of problem that can be dealt with. The recent burst of activity in this field has been 
fueled, in part, by the advent of freely available, efficient implementations of this algorithm 117] . 
There are also interfaces available between the commutative algebra program [T7j and Mathematica 
[18\ [8], with [8] being particularly geared towards physicist's needs. The aim of this talk is to give an 
elementary introduction to the Buchberger algorithm and some of its recent applications. 

In order to give an idea of how one simple algorithm can make so much possible, I will present 
the Buchberger algorithm and then show how it may be applied to physics in 3 elementary examples. 
Firstly, I will describe how it can be used to obtain constraints on the flux parameters in four dimensional 
descriptions of string phenomenological models which are necessary and sufficient for the existence of 
certain types of vacuum [8]. Secondly, I will describe how the Buchberger algorithm can be used to 
simplify the equations describing the vacua of such systems - making problems of finding minima much 
more tractable |8j . Finally, I shall describe how the same simple algorithm can be used to calculate the 
supersymmetric vacuum space geometry of the electroweak sector of the MSSM [lj. 

The remainder of this talk is structured as follows. In the next section, I take a few pages to explain 
the algorithm and the few mathematical concepts that we will require. In the three sections following 
that, I then describe the three examples mentioned above. I shall conclude by making a few final 
comments about the versatility and scaling of the Buchberger algorithm. 

2 A tiny bit of commutative algebra 

Two pages of simple mathematics will suffice to achieve all of the physical goals mentioned in the 
introduction. First of all we define the notion of a polynomial ring. In this paper we will call the fields 
of the physical systems we study (j) 1 and any parameters present, such as flux parameters, a a . The 
polynomial rings C [0*,a Q ] and C [a a ] are then simply the infinite set of all polynomials in the fields 
and parameters, and the infinite set of all polynomials in the parameters respectively. 

Another mathematical concept we will require is that of a monomial ordering. This is simply 
an unambiguous way of stating whether any given monomial is formally bigger than any other given 
monomial. We may denote this in a particular case by saying mi > rri2 where m\,m2 G C [</>', a Q ] are 
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monomials in the fields and parameters. It is important to say what is not meant by this. We are not 
saying that we are taking values of the variables such that the monomial mi is numerically larger than 
the monomial m<i- Rather we are saying that, in our formal ordering, m\ is considered to come before 
m,2. 

For our purposes we will require a special type of monomial ordering called an elimination ordering. 
This means that our formal ordering of monomials has the following property. 



In words this just says that if the largest monomial in P according to our ordering, LM(P), does not 
depend on <p l then P does not depend on the fields at all. The monomial ordering classes all monomials 
with fields in them as being bigger than all of those without such constituents. 

Given this notion of monomial orderings we can now present the one algorithm we will need to use 
- the Buchberger algorithm |15j . The Buchberger algorithm takes as its input a set of polynomials. 
These may be thought of as a system of polynomial equations by the simple expedient of setting all 
of the polynomials to zero. The algorithm returns a new set of polynomials which, when thought of 
as a system of equations in the same way, has the same solution set as the input. The output system, 
however, has several additional useful properties as we will see. 

The Buchberger Algorithm 

1. Start with a set of polynomials, call this set Q. 

2. Choose a monomial ordering with the elimination property described above. 

3. For any pair of polynomials Pi,Pj G Q multiply by monomials and form a difference so as to cancel 
the leading monomials with respect to the monomial ordering: 



4. Perform polynomial long division of S with respect to Q. That is, form h = S — m^P^ where 7713 
is a monomial and P^ G Q such that m^2A{P^) cancels a monomial in S. Repeat until no further 
reduction is possible. Call the result h. 

5. If h = consider the next pair. If h ^ add h to Q and return to step 3. 

The algorithm terminates when all S-polynomials which may be formed reduce to 0. The final set of 
polynomials is called a Grobner basis. 

As mentioned above, the resulting set of polynomials has several nice properties. The feature which 
is often taken as defining is that polynomial long division with respect to this new set of polynomials 
always gives the same answer - it does not matter in which order we divide the polynomials out by. 

For us, however, the important point about our Grobner basis Q is that it has what is called 
the elimination property. Q n C[a a ], the set of all polynomials in Q which depend only upon the 



P G C [<f>\ a a ] , LM(P) G C [a a ] =>PgC [a a ] 



(1) 



S = Pl P I -p 2 Pj s.t. p 1 LM(P I ),p 2 LM(Pj) cancel. 



(2) 
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parameters, gives a complete set of equations on the a a which are necessary and sufficient for the 
existence of a solution to the set of equations we started with. The reason why this is so is actually 
very straightforward. Our elimination ordering says that any monomial with a field in it is greater than 
any monomial only made up of parameters. Looking back at step 3 of the Buchberger algorithm we see 
that we are repeatedly canceling off the leading terms of our polynomials - those containing the fields 
- as much as we can. Thus if it is possible to rearrange our initial equations to get expressions which 
do not depend upon the fields (ft 1 then the Buchberger algorithm will do this for us. Clearly, while we 
have interpreted the a a as parameters and the (ft 1 as fields in the above, as this is what we will require 
for the next section, this was not necessary. The Buchberger algorithm can be used to eliminate any 
unwanted set of variables from a problem, in the manner we have described. 

This completes all of the mathematics we will need for our entire discussion and we may now move 
on to apply what we have learnt. 



3 Constraints 

The first physical question we wish to answer is the following. Given a four dimensional jV = 1 
supergravity describing a flux compactification, what are the constraints on the flux parameters which 
are necessary and sufficient for the existence of a particular kind of vacuum? This question can be 
asked, and answered [8], for any kind of vacuum, but in the interests of concreteness and brevity let us 
restrict ourselves to the simple case of supersymmetric Minkowski vacua. 

Here is the superpotential of a typical system, taken from |19j . It describes a non-geometric com- 
pactification of type IIB string theory. 

W = a - 3a lT + 3a 2 r 2 - a 3 r 3 (3) 
+S(-b + 3&ir - 3b 2 r 2 + 6 3 r 3 ) 
+3U (co + (ci + ci + ci)r - (c 2 + c 2 + c 2 )t 2 - c 3 r 3 ), 

This system has some known constraints on its parameters which are necessary for the existence of a 
permissible vacuum. These come from, for example, tadpole cancellation conditions. 

a &3 - 3ai6 2 + 3a 2 h - a 3 6 = 16 (4) 
aoc 3 + ai(c 2 + c 2 - c 2 ) - a 2 (ci + c x - ci) - a 3 c = 
cob 2 - cih + c\b\ - c 2 b = cqc 2 - c\ + c\c\ - c 2 c = 
ci& 3 - c 2 b 2 + c 2 b 2 - c 3 bi = c 3 ci - c| + c 2 c 2 - cic 3 = 
co^3 - cib 2 + c\b 2 - c 2 b\ = c 3 c - c 2 c\ + c 2 c\ - c\c 2 = 
c\b 2 - c 2 b\ + c 2 bi - c 3 b = c 2 c\ - c\c 2 + c\c 2 - c c 3 = . 

We also have the same constraints with the hats and checks switched around. In this example the fields, 
which we have been calling (ft 1 , are S, t and U and everything else is a "flux" parameter, or an a a in our 
notation. 

In total, the equations which must be satisfied if a supersymmetric Minkowski vacuum is to exist 
are W = 0, dsW = 0, d T W = 0, djjW = and the constraints on the flux parameters given above. 
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To extract a set of constraints solely involving the parameters which are necessary and sufficient for 
the existence of a solution to these equations we simply follow the procedure outlined in the previous 
section. 

We can carry out this calculation trivially in Stringvacua 8 J and, in fact, this example is provided 
for the user in the help system. The result is 

= C\ = C2 = C\ = C2 = Cl = C2 = Co = C 3 (5) 

= 16 + 0360 — 3a 2 bi + 3ai62 — ao&3 

= 16a§6§ - 96a 2 a 3 6 6i - 288a|6f + 432aia 3 6? + 54a^f - 81aia 2 a 3 6? + 27a a^6? + 432aia 3 6 6 2 

OOOO o 00 

— 27a 2 a 3 6 62 + 48aia 3 6 62 — 288aoa 3 6i62 — lSaia 2 a 3 bobib 2 — 45aoa 3 6o^i^2 — 54aia 2 b 1 b2 

O O O O Q 

+81a 1 a 3 6 1 62 — 27aoa 2 a 3 b 1 b 2 + 54aoa2a3&o^2 + 27aoaia 3 &i b 2 — 27a a 3 6 2 — 288ai 02^0^3 

Q O - o o 

— 32aoa 3 6 6 3 + 27a 2 b 6 3 — 45a 1 a 2 a 3 6 6 3 + 432a a2&ib 3 — 27aia 2 b bib 3 + 54a 1 a 3 6 &i6 3 

o o o 

+48a a 2 a 3 6 6i6 3 + 18a a 2 6 1 6 3 - 81a aia 3 6 1 6 3 - 144a ai6 2 ^3 + 27a 1 a 2 b b 2 b 3 - 54a a 2 b b 2 b 3 
-51a aia 3 b b 2 b3 + 27aoaia 2 6i6 2 6 3 + Aba^a 3 bib 2 b 3 - 27a alblb 3 + 27a^a 2 b 2 b 3 + 16aQ6| 
— 27af 60^3 + 45aoaia26o^i + 27aQa\b\b\ — 48aga26i6 3 + 2>a^a\b 2 b\ 

The reader will note that the result is a somewhat lengthy set of equations. In principle one has 
to find quantized solutions to these expressions, an obviously intractable Diophantine problem, and 
therefore it might be asked why this result is of any use. In fact, knowledge of such constraints on the 
flux parameters is hugely useful for a number of reasons. 

• Firstly, we note that, while the full result of this process is often complex, some of the constraints 
can give us simple information about the system. In the current case, for example it can be seen 
that c 2 = is required for the existence of a supersymmetric Minkowski vacuum. 

• Secondly, if one is scanning over a range of flux parameters and trying to numerically solve the 
equations to find vacua one can speed up ones analysis by first substituting any given set of flux 
parameters into the constraints we have obtained. If the constraints are not satisfied then vacua 
do not exist and there is no point in searching numerically for a solution. This turns what would 
be a time consuming numerical process giving inconclusive results (no solution was found) into a 
quick analytic conclusion (no solution exists). 

• Lastly, knowledge of such constraints can greatly speed up algebraic approaches to finding vacua 
such as those outlined in [8|. 

4 Simplifying equations for vacua 

Another use for the mathematics we learnt in section [2] are the so called "splitting tools" used in work 
such as [8]. The physical idea here is simple. Consider trying to solve the equations dV/dcj) 1 = to 
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find the vacua, including those which spontaneously break supersymmetry, of some supergravity theory. 
These equations are often extremely complicated. One way of viewing why this is so is that the equations 
for the turning points of the potential contain a lot of information. They describe not only the isolated 
minima of the potential which are of interest, but also lines of maxima, saddle points of various sorts 
and so forth. A useful tool to have, therefore, would be an algorithm that takes such a system as an 
input and returns a whole series of separate sets of equations, each individually describing fewer turning 
points. Since each separate equation system would then contain less information one might expect them 
to be easier to solve. It would be beneficial to choose a division of these equations which has physical 
interest. The choice we will make here, and which programs like Stringvacua implement [8], is to split 
up the equations for the turning points according to how they break supersymmetry - that is according 
to which F-terms vanish when evaluated on those loci. 

The ability that packages such as Stringvacua have to split up equations in this manner is based upon 
the following splitting tool (see [21] for a nice set of more detailed notes on this kind of mathematical 
technique). Say that one of the F-terms of our theory is called F. Then we can split the equations 
describing turning points of the potential into two pieces. 



The first of these expressions is a set of equations which is easier to solve, in general, than dV/d^ 1 = 
alone. We can use the F-term to simplify the equations for the turning points of the potential. On the 
other hand, expression ([7]) is not even a set of equations - it contains an inequality. We can convert ([7]) 
into a system purely involving equalities by making use of the mathematics we learned in section (2) 

Consider the following set of equations, including a dummy variable t. 



The second equation in ([8]) has a solution if and only if F ^ 0, simply t = 1/F. If F = the equation 
reduces to — 1 = which clearly has no solutions. The equations ([8]), then, have a solution whenever the 
set of equalities and inequalities (J7J do. Unfortunately they also depend upon one extra, and unwanted, 
variable - t. This is not a problem as we already know how to remove unwanted variables from our 
equations. We can simply eliminate them, as we did the fields in section [2j This will leave us with a 
necessary and sufficient set of equations in <j> 1 and a a for a solution to ([8]) and thus to dTJ). 

So we can split the equations for the turning points of our potential into two simpler systems. One 
describes the turning points of V for which F = and the other those for which F / 0. We can of 
course perform such a splitting many times - once for each F-term! In fact, using additional techniques 
from algorithmic algebraic geometry \20\ [8]. which are essentially based upon the same trick, one can 
go much further. One can split the equations for the turning points up into component parts gaining 
one set of equations for every separate locus. Because we know which F-terms are non-zero on each of 
them these are classified according to how they break supersymmetry. The researcher interested in a 
certain type of breaking can therefore select the equations describing the vacua of interest and throw 
everything else away. 



dV/d(j) 1 = , F = 
dV/dfi = , F ^ 



(6) 
(7) 



dV/d<f> % = , Ft - 1 = 



(8) 
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The above process of splitting up the equations for the vacua of a system can be very simply carried 
out in Stringvacua. Numerous examples can be found in the Mathematica help files which come with 
the package [5]. Here, let us consider the example of M-theory compactified on the coset 777^7777^- 
The Kahler and superpotential for this coset, which has SU(3) structure, has been presented in [22J. 

k = -4io g (- i ([/-c7))-iog(-i(r 1 -T 1 )(r 2 -r 2 )(r3-T3)) (9) 

W = -= [4U(T X + T 2 + T 3 ) + 2T 2 T 3 - T X T 3 - T X T 2 + 200] 

Even this, relatively simple, model results in a potential of considerable size. Defining Tj = — it{ + r, 
and U = —ix + y we find, 

v = K^7 7 7 4 ( 40000 + l l T l - 400tit 2 - 4t|nr 2 + 4^r| + rfrf - 400tit 3 + 800r 2 r 3 + 

2r 1 2 r 2 r 3 - 4r 1 r 2 2 r 3 + t 2 t 3 2 - 4nr 2 r 3 2 + 4t 2 2 t 3 2 - 24t 2 t 3 x 2 + 4t|x 2 - 24*i (i 2 + t 3 )x 2 (10) 
+4r 2 x 2 + 8rir 2 x 2 + 4r|x 2 + 8tit 3 x 2 + 8r 2 r 3 x 2 + 4r|x 2 + 1600riy - 8tlny 
+1600r 2 y + l6t 2 3 T 2 y - 8rfT 2 y - Sr^y + 1600r 3 y - 8rfT 3 y + 16r|r 3 y - 8rir 3 2 y 
+ 16r 2 r 2 y + 16t|y 2 + 16r 2 y 2 + 32r 1 r 2 y 2 + 16r|y 2 + 32r 1 r 3 y 2 + 32r 2 r 3 y 2 + 16r|y 2 

+ *i + r 2 + 2r 2 r 3 + rf + 4x 2 - 8r 2 y - 8r 3 y + 16y 2 ) + ^(4t 2 + rf - 4ri(r 3 + 2y) 
+4(r 3 2 + x 2 + 4r 3 y + 4y 2 )) . 

To find the turning points of this potential we naively need to take eight different derivatives of (|10p 
and solve the resulting set of inter-coupled equations in eight variables. This is clearly prohibitively 
difficult. Using the techniques described in this section, however, Stringvacua, can separate off parts of 
the vacuum space for us with ease. Consider, for example, the vacua which are isolated in field space 
and for which the real parts of all of the F-terms are non-zero, with the imaginary parts vanishing. To 
find these, the package tells us, we need only solve the equations, 

9x 2 - 500 = , 5ii - 2x = , t 2 - x = , t 3 - x = , n = r 2 = r 3 = y = . (11) 

Because they only describe a small subset of all of the turning points of the full potential these equations 
are extremely simple in form and may be trivially solved. For this particular example the physically 
acceptable turning point that results is a saddle - something which can be readily ascertained once its 
location has been discovered. 



5 Geometry of vacuum spaces 

As a final example of what we can do with the simple techniques introduced in section [2] we will show 
how to calculate the vacuum space of a globally supersymmetric gauge theory. It is a well known result 
(see [23] and references therein) that the supersymmetric vacuum space of a such a theory, with gauge 
group G, can be described as the space of holomorphic gauge invariant operators (GIO's) built out of 
F-flat field configurations. What does this space look like? Consider a space, the coordinates of which 
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Type 


Explicit Sum 


Index 


Number 


LH 




i = 1,2,3 


3 


HH 






1 


LLe 




i,j = l,2,3;k = l,...,j-l 


9 


LHe 




i,j = 1,2,3 


9 


V 


v i 


i = 1,2,3 


1 



Table 1: The set of elementary gauge invariant operators for the electroweak sector of the 
MSSM. 

are identified with the GIO's of the theory. If there were no relations amongst the gauge invariant 
operators then this space would be the vacuum space. However there frequently are relations because 
of the way in which the GIO's are built out of the fields. For example, if we have three gauge invariant 
operators S 1 , S 2 and S 3 which are built out of the fields as S 1 = (0 1 ) 2 , S 2 = ((f) 2 ) 2 , S 3 = (p 1 ^ 2 then we 
have the relation S^^S 2 = (S 3 ) 2 . If we take these GIO's to be built out of the F-flat field configurations 
then there will be still further relations among them. The vacuum space of the theory is the subspace 
defined by the solutions of these equations describing relations amongst the gauge invariant operators, 
once F-flatness has been taken into account. 

How can we calculate such a thing? The holomorphic gauge invariant operators of a globally 
supersymmetric gauge theory are given in terms of the fields. 

S 1 = /V) (12) 

Here S 1 are our GIO's and the f 1 are the functions of the fields that define them. Let us write the 
F-terms of the theory as F' 1 . Consider the following set of equations. 

F i = , S 1 - f 1 (</>*) = (13) 

These equations have solutions whenever the S 1 are given by functions of the fields in the correct 
way and when those field configurations which are being used are F-flat. However, according to the 
proceeding discussion, we wish to simply have equations in terms of the GIO's to describe our vacuum 
space. As in previous sections, we can eliminate the unwanted variables in our problem, in this case the 
fields (p' 1 , using the algorithm of section [5] to obtain the equations describing the vacuum space. 

As a simple example, let us take the electroweak sector of the MSSM [lj (with right handed neutri- 
nos). Given the field content of the left handed leptons, L l a , the right handed leptons, e l and i/, and 
the two Higgs, H and H, one can build the elementary GIO's given in table [TJ The indices i,j run over 
the 3 flavours and the indices a, (3 label the fundamental of SU(2). 

To compute the F-terms we require the superpotential. Let us take the most general renormalizable 
form which is compatible with the symmetries of the theory and R-parity. 

^minimal = C° £ H a H^ a/3 + ^ C 3 ^ £ I^H^ + ]T C^V + £ C%V l £ VJl^ . (14) 

a,j3 i,j a,/3 i,j i a,/3 
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Here e is the invariant tensor of SU(2) and C°, Cfp C^-, and, Cfj are constant coefficients. 

We now just follow the procedure outlined at the start of this section. We calculate the F-terms by 
taking derivatives of the superpotential, we label the gauge invariant operators Si to £23, we form the 
equations (|13p and then simply run the elimination algorithm given in section [2j 

The result is, upon simplification, given by six quadratic equations in 6 variables. It is a simple 
description of an affine version of a famous algebraic variety - the Veronese surface pp. What can be 
done with such a result? The first observation we can make is that this vacuum space is not a Calabi- 
Yau. This means, for example, that one can say definitively that it is not possible to engineer this 
theory by placing a single D3 brane on a singularity in a Calabi-Yau manifold, without having to get 
into any details of model building. 

Secondly one can study such vacuum spaces in the hope of finding hints at the structure of the 
theory's higher energy origins. In the case we have studied in this section, for example, we can "pro- 
jectivize" (pretend the GIO's are homogeneous coordinates on projective space rather than flat space 
coordinates) and study the Hodge diamond of the result. The structure of supersymmetric field theory 
tells us that this Hodge diamond should depend on 4 arbitrary integers, but there is nothing at low 
energies which prevents us from building theories with any such integers we like. Interestingly, in the 
case of electroweak theory, these integers are all zero or one. 

h°>° 1 
/i ' 1 /i ' 1 

h P,q = h 0,2 h l,l h 0,2 _ 1 . (15) 

/i ' 1 /i ' 1 

h°>° 1 

Whether this structure is indeed a hint of some high energy antecedent or just a reflection of the 
simplicity of the theory is debatable. This example does, however, demonstrate the idea of searching 
for such evidence of new physics in vacuum space structure. We should also add here that similar 
techniques can be used to show that the vacuum space of SQCD is a Calabi-Yau [5]. 



6 Final Comments 

To conclude we shall make several points - one of which is a note of caution, with the rest being more 
optimistic. The first point which we shall make is that we should be careful lest the above discussion 
makes the algorithm we have been describing sound like an all-powerful tool. There is, as ever, a catch. 
In this case it is the way the algorithm scales with the complexity of the problem. A "worst case" upper 
bound for the degree of the polynomials in a reduced Grobner basis can be found in [23]. If d is the 
largest degree found in your original set of equations then this bound is, 

2 {j + d ) ' (16) 

where n is the number of variables. This worst case bound is therefore scaling doubly exponentially 
in the number of degrees of freedom. These very high degree polynomials are an indication that the 

8 



problem is becoming very complex and thus computationally intensive. Despite this, physically useful 
cases can be analysed using this algorithm quickly, as demonstrated in this talk and in the references. 
This scaling does mean that one is not likely to gain much by putting one's problem on a much faster 
computer. One good point about (fT6|) is that if you can find a way, using physical insight, to simplify 
the problem under study, then what you can achieve may improve doubly exponentially. Such a piece 
of physical insight was one of the keystones of the application of these methods to finding flux vacua 

We finish by commenting that the methods of computational commutative algebra which we have 
discussed here are extremely versatile. We have been able to perform three very different tasks simply 
utilizing one algorithm in a very simple manner. These methods are of great utility in problems taken 
from the literature and their implementation in a user friendly way in Stringvacua means that they 
may be tried out on any given problem with very little expenditure of time and effort by the researcher. 
Many more techniques from the field of algorithmic commutative algebra could be applied to physical 
systems than those described here, or indeed in the physics literature. We can therefore expect that 
this subject will only increase in importance in the future. 
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